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' Abstract 

\^ • A generalized theory of two-dimensional isotropic turbulence is developed based on con- 



formal symmetry. A number of minimal models of conformal turbulence are solved under an 



' extended constraint including both the enstrophy cascade by Kraichnan and the discontinuity 



of vorticity by Saffman. There are an infinite number of solutions which fall into two different 
0^ ■ categories. An explicit relation is derived in one of the categories between the central charge 

of Virasoro algebra, the lowest anomalous dimension and the power of the energy spectrum. 
^ • Some statistical properties such as energy spectrum, skewness, flatness and Casimir invariants 

. are predicted and compared with numerical simulation by the pseudospectral method. 
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1 Introduction 

The problem of the energy spectrum has still been controversial in two-dimensional turbulence. In 



three dimensions, dimensional analysis by Kolmogorov |1£] shows that the energy spectrum in the 



inertial subrange should be scaled as E{k) = Ce^^^k where e is the energy dissipation rate. 



Kraichnan |17], Leith |2C] and Batchelor ||l| applied Kolmogorov theory to the two-dimensional 
turbulence and obtained that a power of the energy spectrum should be —3 according to existence 
of the enstrophy cascade in the inviscid limit. On the other hand, Saffman derived the 
power —4 of the energy spectrum by assuming the jump of constant- vorticity layers. A possible 



explanation of these difference is proposed by Gilbert |12| in which spiral structures of vortical 



region play an important role and may fill the gap between these two spectra. 
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In many of numerical simulations [16|[P|2^ of decaying two-dimensional turbulence with the 
initial gaussian spectrum, the power increases from —4 to —3, contours of vorticity are stretched 
and an active enstrophy cascade is observed up to t ~ 10. Then coherent vortices emerge with 
spiral tails and the power of the energy spectrum decreases [^]. Moreover, a very long-time 
simulation [^ ] shows that in the final stage the flow seems to settle down to the universal "sinh" 
state with the maximum entropy, which is predicted by a statistical argument of the assembly of 
point vortices [^. Analytical closure theories, such as direct interaction approximation (DIA) 
and test field model have been also applied to evaluate quantitatively the energy spectrum with its 
constant and the enstrophy transfer |18| in a reasonable comparison with numerical results. 
Kraichnan [^] argues that the logarithmic correction can escape from the divergence of the 
enstrophy integral. 



Recently, Polyakov |26| |27] has proposed that one may use the conformal field theory (CFT) 
to analyze two dimensional turbulence. Since there appears power law in the correlation functions 
of velocity fields, for example, the existence of the conformal invariance seems quite natural. The 
merit to use CFT is that there are well-classified models (minimal models) where there appears 
only finite number of independent operators. Usually, when we try to calculate n point functions, 
we need the knowledge of n + 1 point functions because of the nonlinearity of the Navier-Stokes 
equation. Since we can not usually cope with infinite number of degree of freedom, we need 
to restrict them to the physically relevant ones from outside, for example, quasi-normal and 
DIA theories. If the turbulence is described by CFT minimal models, such restriction comes 
from mathematical consistency and there are no arbitrariness. In the original paper, Polyakov 
proposed a minimal model (21,2) as the candidate which may describe turbulence. This model 
predicts quite a reasonable value for the exponent. However, it was indicated later by several 



groups 1 22 1 [ 21 1 [Q] 1 11] that Polyakov's condition on the CFT is not satisfied by unique model, but 
in fact there are infinite number of candidates. Also, as Polyakov himself discovered |27|, there 
may be some ambiguities in the conditions for CFT if one considers the non-vanishing of one 
point functions in the system with boundaries. Even if one takes the modified conditions, the 
situation does not change and there are again a lot of CFT candidates. Q In any case, if one can 
not find other appropriate conditions, one may not tell which model is more preferable. One can 
not even say whether there are several universality classes. 

Due to appearance of infinite solutions, we can not predict uniquely the exponent of energy 
spectrum. In such situation, it is difficult to convince ourselves whether prediction from CFT 
is useful as long as we consider only one observable. One way to proceed in this situation is 
to consider the global (or geometrical) aspect of the conformal turbulence. |5| In this paper, we 
would like consider another direction. In CFT models, there appear a lot of primary fields beside 
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the original velocity field. It is quite important consider the physical meaning of these fields. 
Some of them can be obtained in the short distance expansion of the velocity field. Thus we can 
relate it to the the dependence of three or four point functions on the Reynolds number which 
is known as the "skewness" or "flatness". If the assumption of the local equilibrium is correct, 
they should be some universal constants. However, in CFT, we get nontrivial dependence on 
Reynolds number due to the anomalous dimension. In three-dimensional turbulence, there is a 
tendency that both skewness and flatness of velocity derivatives increases as Reynolds number 
becomes large [^]. In many of CFT solutions, however, the behavior is the opposite to these 
three-dimensional observations. 

Interestingly, as we prove in this paper, these exponents are not independent with each other 
but can be uniquely determined from the exponent of energy spectrum. This property is tightly 
related to the finiteness of minimal models. Intuitively, in order to calculate the higher correlation 
functions, the most dominant contribution comes from the lowest dimension operator in the OPE. 
However, in the minimal models, there are always lower bound in the dimensions of primary fields. 
We surveyed more than 1600 candidates of CFT solutions. In 99% of them, the lowest dimension 
is attained in the first OPE between ip x ij:. The other 1% behaves quite differently but we show 
that they can be determined exactly. 

As another peculiar character, two-dimensional inviscid flows possess infinite number of con- 
served quantities < > which are topological invariants of area-preserving diffeomorphism, 
called Casimir. It is commonly believed that these conservation laws may affect the character of 
the two-dimensional turbulence, i.e. ergordicity, energy spectrum and so on, although the effects 
are not fully understood. A finite-mode system of two-dimensional Euler flow preserving such 



conservation laws is derived by Zeitlin[32|. However, the effect of this conservation laws on the 
behavior of turbulent states is still unresolved. Using CFT theories, Falkovich and Hanany[P 
showed that the dependence of the dissipation rate of < > on the viscosity v can be predicted. 
We extend their results and compare the prediction with the numerical results. 

This paper is organized as follows. In section 2 we review the classical theories of two- 
dimensional turbulence. It is shown that the extended theory, which contains both the Kraichnan- 
Leith-Batchelor's and Saffman's theories, is useful in order to understand the conformal turbu- 
lence. In section 3 Polyakov's conformal turbulence is reviewed with a brief sketch of conformal 
field theory. The original CFT models are generalized so that the imposed constraint includes not 
only the Kraichnan's constant transfer of enstrophy but also the non-cascade theory by Saffman. 
In section 4 solutions of this generalized CFT models of two-dimensional turbulence are shown. 
Distributions of the energy spectra are revealed, A relation between the central charge in CFT 
and the power of the energy spectrum is described. In section 5 higher order statistics including 
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skewness and flatness factors are derived. Scaling of Casimir invariant < cj" > and its decay 
rate are discussed in section 6. These predictions are compared with the numerical simulations 
in section 7. Conclusions and discussions are summarized in section 8. Throughout this paper, 
attention is paid on removal of the gap between the field theory and fluid communities. 



2 Theories of two-dimensional turbulence 

We shall review standard theories of two-dimensional turbulence in order to understand the 
conformal theory of it. In what follows we mainly keep in mind the decaying turbulence, i.e. the 
system with no external forcing. The equations governing the two-dimensional incompressible 
viscous flow are 

dv , , 

— • V)v = -Vp + uAv, (1) 
V-v = 0. (2) 

The vorticity equation is given by 

^ + J{i;,uj)=uAuj, (3) 
where ■0 is the streamfunction, uj = —Aip is the vorticity, and the Jacobian J is defined by 

' dy dx dx dy ' 

Prom (H) and (|3|), the evolution of total energy E = ^ < Iv]"^ > and square- vorticity, called 
enstrophy Q = ^ < >, is given by 

f = -2.^^-e. (5) 

§ = -2.P.-,. (6) 

Because the enstrophy decreases monotonically, the energy dissipation rate e vanishes in the 
inviscid limit u ^ 0. Thus the stationary energy cascade is irrelevant in decaying two-dimensional 
turbulence, although in the forced system, the energy is transferred to the low wavenumber side, 
called the inverse cascade. 

In Kraichnan-Leith-Batchelor theory, the turbulence is characterized by the cascade of en- 
strophy. The dimensional analysis shows the energy spectrum is expressed as 

E{k) = Cri^l^k-^, (7) 
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in the inertial subrange L^^ <^ k <^ 6^^ where rj is the enstrophy dissipation rate, C is a universal 
constant, L denotes the large scale length which may be considered as the size of the largest eddy 
of coherent structures or the boundary condition, and 5 implies the dissipation scale given by 

The total energy estimated by the integral in the inertial range leads to 

E = J ^ E{k)dk = -Crj^/^L^. (9) 
On the other hand L may be evaluated by E and as 

L-i^^f. (10) 



Prom (I) and (|^) we may obtain the Reynolds number as 



Re^^^ t) ■ (11) 



In this theory rj and v are considered as independent and it is allowed that rj remains finite in 
the inviscid limit. 

Saffman considered that, in the context of decaying two-dimensional turbulence, the discon- 
tinuity of the vorticity is the asymptotic feature by extending the argument of one-dimensional 
turbulence, i.e. Burgers equation. The energy spectrum may be given by 

Eik) = -^k-^, L-^<^k<^5-\ (12) 



and the dissipation length is given by 



The enstrophy dissipation rate is estimated as 



^2/3^7/6 



Using LH} = E and E remains constant in the inviscid limit, (14) is rewritten as 

^2/35^11/6 

Thus rj depends on v as v"^/^ . By Saffman's picture there is no enstrophy cascade in the inviscid 
limit. Reynolds number is given by 

Re=U). (10) 
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Comparing the energy spectra and the dissipation lengths of the two theories, it seems quite 
useful to develop a following extended theory. That is, the energy spectrum in the inertial range 
is given by 

^(fc) = L-^<^k<^6~\ (17) 

and the dissipation scale is expressed as 

6 = l(^^Y, (18) 

where parameters a and /? are introduced in order to bridge over the two standard theories. This 
generalization is also important when we consider CFT turbulence described in the following 
section. 



From (18), the Reynolds number is estimated as 

The enstrophy dissipation rate rj is estimated as 

V = 2u j^^^ k'E{k)dk^^-^. (20) 
Substituting (|T^) into (pO[ ) the dependence of on 17, zv, L is expressed as 

r? « f]i+^/2j.i-7L-2+27^ (21) 



where 7 = /3(a + 5). Using ( pUj ) rj is expressed by (7, z^, and E as 

?L^v / E 

Moreover, the restrictions 

(a,/3) = (-3,^), (7 = 1), (a,/3) = (-4,i), (7 = i), (23) 

may be imposed if this extension should meet both the enstrophy cascade theory and Saffman's 
spectrum. 

A more physical argument, which reconcile the two theories similarly to our idea, is made 



by Gilbert |12|, who showed that the energy spectra with powers between —3 and —4 can be 
derived by considering that the isolated strong vortex cores wind up the discontinuous lines of 
the vorticity. This spiral model gives the energy spectrum in the inertial subrange as 

E{k) « (rt)2/(^+i) k-^+{s-i)/(s+i) ^ (24) 

where T denotes the strength of the vortex and the parameter s specifies the distribution of the 
vorticity in the core. The value of s is allowed as 1 < s < 00, s = 1 implies the Saffman's power 
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—4, s = oo gives the KLB values —3 and the index of the reahstic vortex is approximately 2 ( 
the value of the point vortex ), which gives the power —11/3. 

In the case of decaying turbulence, the temporal behavior of 17 can be obtained by integrating 
(|2^) assuming that E is nearly constant: 
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o(t) = [fi(or^+^/2 + (i-|)(^) ""^]-^ (25) 

Thus the enstrophy decays at large t as 

17«(|)^t^, (26) 
and the enstrophy dissipation rate becomes 

7(1-7) 4_ 

ry^Q-^t^. (27) 
These are consistent with ^ t^'^,T] ~ t^'^ of Batchelor if 7 = 1 and Q ^ t^^'^,r] ~ 

t-lV5 of 

Saffman if 7 = 1/3. 

We are also interested in the behavior of the infinite number of inviscid invariants, i.e., ar- 
bitrary powers of vorticity, in the limit u ^ 0. The decay rate of i7„ = ^ < cj" > is given 

by 

di}„ n(n — l)^ „ 9,„ ^9 , ^ 

Vn ^ = ^ ^ ' < u:^-\VLof > . (28) 

It is difficult to get the dependence of r]n on Q and t], but if we assume a random phase 
distribution of the Fourier coefficients of w, we get the dependence as 

rj^ = 2t-in(n - l)(n - 3)!!i/0t-ip ^ n^-^rj, (29) 
for an even n and t/„ = r2„ = for an odd n. From (|2^) , (^) and (|2^ , rjn decays as 

(l-7)["-2+7] 



, , n-\-2 — 'y 

r?nOC - t^. (30) 



The pair {a, 13) = (—3,1/2) gives 
and (a,/3) = (—4, 1/3) leads to 



rjn ~ (i?/z,)(6--10)/15^-l~3n/5_ 

Thus the dissipation rate of is independent of viscosity by the enstrophy cascade theory and 
increases as v tends to zero by Saffman's picture in the random phase approximation. 
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3 Conformal turbulence 



3.1 Short summary of CFT 

In order to make smooth communication between fluid dynamics and elementary particle com- 
munity, we would like give a short sketch of some basics of CFT. The material is quite standard. 
See the original paper to understand the detail[||. Also there are many nice review articles, for 
example, see ||T5|| . 

In quantum field theory, symmetry generally plays an important role. On the one hand 
symmetry implies existence of conservation law and resulting in effective reduction of the degrees 
of freedom. On the other hand, symmetry is a useful tool to select out nicely behaving class of 
field theories among vast amount of all the possible field theories. 

The most successful example of such is the class of two dimensional conformal field theories 
(CFT). The field theory in the class is defined in two dimension, and has the conformal invariance 
as its symmetry. This means, if we use a complex variable z = xi + ix2, dynamics of theory is 
invariant under an complex analytic transformation z' = f{z). Polyakov et al. performed a 
systematic study of CFT and showed that the conformal symmetry is really powerful to classify 
possible CFTs and solve them. 

Physically, CFT concerns the critical point of the second order phase transition of statistical 
systems. On the critical point, we are faced with the situation that the correlation length of 
the operators in the system grows up to infinity. In this limit, the system gains the conformal 
invariance. Then the two dimensional system is described by CFTs. 

Especially in two dimensions, it is known that the conformal symmetry becomes infinite 
dimensional and generated by the operators L„ and Z„ (n = 0, ±1, ±2, • • •). They form a Lie 
algebra: 

[Ln, Lm] = (n- m)Ln+m + ]^("'^ ~ n)5n+mfl, (31) 

which is called Virasoro algebra. The last term in the RHS comes from quantum anomaly. 

Important fact is that we have a detailed understanding of its representations theory, and he 
knowledge is really helpful in classifying all the possible CFTs. Representation theory is quite 
different whether c > 1 or c < 1. In the former case, the classification of the possible theories 
is not completely done with Virasoro algebra alone, we need further informations such as the 
extended symmetry to classify and solve the system. Conversely, in the latter case, the Virasoro 
algebra is strong enough to completely classify all the possible theories and to solve the system 
completely. In fact, it was shown that the central charge is restricted to take only the discrete 
values, 

c=l-5<^^, (32) 

pq 
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with {p, q) mutually coprime integers. These models are so-called minimal models. 

In CFT, we have two kinds of fields called primary fields and secondary fields. The former is 
covariant under the conformal transform, 

(Paa{z^z) ^ [-g^ j [-g^j (I>aa{w,w), (33) 

while the latter is not. For each primary field, there are infinite number of secondary fields. Those 
can be obtained from the primary field through the action of the Virasoro generators. Derivative 
of the primary field is a typical secondary field. 

One of the most important properties of the minimal model is that the number of primary 
fields is finite. For the minimal model characterized by p, q, the conformal dimensions A (A) of 
primary fields are restricted to take the following discrete values, 

Ars = ~ l<r<p-l, l<s<q-l, ps>qr. (34) 

4pq 

The final restriction ps > qr comes from the fact that we can identify primary field with label 
(r, s) and {p — r,q — s). This set of primary fields is called the Kac table. 

It may be intuitive to compare the representation of the Virasoro algebra with that of su{2) 
algebra]^ The primary field described above corresponds to the highest weight state in the su{2) 
representation. The states generated from highest weight state by the action of lowering operator 
should be compared to the secondary field. Finally, corresponding to the composition rule of the 
angular momenta, two primary fields are composed to give several other primary fields. Namely 
the primary fields form an algebra by the composition operation. The composition rule is called 
fusion algebra or operator product expansion, OPE. 

The rule can be deduced from the representation theory just like the computation of the 
Clebsh-Gordan coefficients of su{2) representation. For the Virasoro algebra, the OPE is explicitly 
given by, 

miTi(r' +r" ~l,2p—r' —r" — 1) min{s' +s" ~l,2q—s' —s" — 1) 

(j)r's' X (t>r"s" ~ X! X! t'^''*] ' (^^) 

r=\r'-r"\+l s=\s'-s"\+l 

where (/>r,s denotes the primary field whose conformal dimension is Ars- The symbol [ ] denotes 

the conformal family (the primary field and its associated secondary fields). If we regard (|35| ) 

as the short distance expansion, the contribution from the secondary fields are higher order in 

power of zz compared to those coming from the primary field, so they are irrelevant. 

Finally, let us see some examples of CFT. The minimal models are not necessarily unitary. 

In other words, negative norm state appears in some situations. Only ifp = q+ l = m+ l with 

^ There are no object in su{2) algebra which correspond to the central charge of the Virasoro algebra. This 
term originates from the infinite dimensionality of the algebra. 
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m = 2,3,4..., the theory is unitary. Many of the critical points of the known 2-dim statistical 
systems are described by the models in this series. For example, the Ising model (c = 1/2, i.e. 
m = 3) the tri-critical Ising model (c = 7/10, i.e. m = 4), the three state Potts model (c = 4/5, 
i.e. m = 5), and so on. For these models, one can find the physical interpretation for the primary 
fields. For example, for Ising model, we have three independent primary operators, 

An=0 A2i = l/16 A3i = l/2. 

These operators may be identified with, identity operator, energy density and the spin field 
respectively. 

On the other hand, the nonunitary models are by no means physical, but contains some very 
interesting models, such as Yang-Lee edge singularity (c = —22/5, i.e. {p,q) = (5,2)). They 
are characterized by the appearance of the negative dimension operators. For the Yang-Lcc edge 
singularity case, we have A2,i = —1/5. In the unitary case, the lowest dimension operator in the 
system is always the vacuum, which is represented by the primary field with Ai^i = 0. For the 
non-unitary case, this 'vacuum' is no longer the lowest dimension operator but rather we have a 
primary field with dimension. 

This operator may be regarded as the 'true vacuum' of the system. 

Let us take a nonunitary model {p,q) = (21,2) as an example. Field content consists of ten 
fields: (pi^i, 02,1, • • • , 010,1- Since the second index s is always unity in this case, we will omit it 
hereafter. OPE is, for example 

04 X 04 ~ [07] + [05] + [03] + [0l] • (37) 

For the two point function, this fusion rule implies that 

< 04(^2)04(0) > = (zz)'^^-^^4 < 07(0) >+••• + {zzf^'-^^^ < 01 (0) > 

+ higher order terms w.r.t. zz. (38) 

Note that z, z dependence of each term in RHS is determined by the accordance of the length 
dimension in both sides. 

3.2 Conformal Turbulence 

3.2.1 Polyakov's original definition 

In CFT approach to two dimensional turbulence, the fluid circulation function tp is supposed to 
be described by a primary field of a minimal model of CFT. ip is related to the velocity field 
through Va = ea/sOpip. Polyakov's conditions for CFT to describe turbulence are following. 
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1. Hopf equation should be satisfied, 



< e^pdai^d'df^i; • • • >= 0. (39) 

2. Constant flow of enstrophy in the momentum. 
Assuming operator product expansion (OPE), 

m [V] ~ [0] + • • • , (40) 

where (p is the smallest dimension primary field in OPE, the nonlinear term in the Navier-Stokes 
equation (5^) is equivalent to 



(L_2i-i - L^2L\)4>, (41) 



where a is a viscous cutoff length. Vanishing of this quantity then produces either 

2A^ - 4A^ > 0, (42) 

or L_20 ~ -^-10- The second equation (which is equivalent to that (f) is (2,1) type primary field) is 
not compatible with the conditions on enstrophy and Polyakov has picked up the first condition. 

The constant enstrophy flow constraint is more subtle. The time derivative of the enstrophy 
gives 

^ =< >, (43) 

with uj = d'^ip and (j)^"^^ = (L_2-L^]^ — L_2-L?_x)</'- Rewriting the right hand side of ( |4^ ) as the 
integral in the momentum space, the constant enstrophy condition is equivalent to 

< tj(x)(/.(2)(0) >~ 3;°. (44) 

Naive dimensional counting shows that it is equivalent to + A^ + 3 = 0. However, we have 
to be careful to the contribution of the one point function, which is non-vanishing in the general 
topology of two dimensional space. The dependence of one point function on the characteristic 
length of the system is 

<V;>=L-2A^, (45) 

for any primary field if). For unitary theory where dimension of any primary field is positive, one 
point function is not relevant in the large L limit. However, since we are discussing on the system 
where dissipation takes place, we are dealing with non-unitary theory where negative dimensional 
operator appears. In this situation, the lowest dimension operator becomes most relevant one. If 
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we denote x the lowest dimension primary field in OPE [i/j] [i;^], then the correct condition for 
the constant enstrophy flow is 

+ A0 - + 3 = 0. (46) 

The solutions of CFT that satisfy (^) and ( |4^ ) was discussed by Chung et. al.Q and they 
found several hundreds candidates. 

Once we know that there are solutions, the next step is how to extract physical exponents out 
of them. The most important one is discussed by Polyakov, i.e. the exponent of energy spectrum, 
which may be described by, 

r 1 2 

E= dkE{k) = -J2 < ^^^^ > • (47) 



Prom contribution of the one-point function (|45|), one finds that 

E{k)r^k", with a = 4A^ - 2A<^ + 1. (48) 
If the one-point correlation integral vanishes, the energy spectrum becomes 

E{k) ~ with a = 4A^ + 1. (49) 

These expressions are Fourier transform of two-point functions, 

< V'(z)V'(0) >~ (50) 
for correlation function on sphere Q and 

< i^{z)i^{0) >~ L-2A0|^|-4A^+2A^^ ^^^^ 

for correlation function on general topology Riemann surface including torus. This kind of corre- 
lation functions was discussed previously by J. Cardy Q which is the refinement of the discussion 
given by Pisher Q. The correlation function of unitary model (in particular, Ising model) on 
the torus was obtained explicitly by Di Francesco et. al.pH]. In the turbulence context, some of 



the physical implication was discussed by Polyakov|27] and Chung et.al. |5[ and [Q. The most 
important observation is that we can not neglect the contribution of the one point function. This 
is because the dimension of (p is negative and it becomes dominant in the large L limit in (|5ll). 

3.2.2 Generalized conformal turbulence 

Let us introduce the other possibilities that can be imposed on the conformal field theory. As 
we saw in section 2, there are many variants of the two-dimensional turbulence models where 
the enstrophy flow may depends on the viscosity. Also, the condition will be changed if one 
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have vanishing one point function which was originaUy considered by Polyakov. To get the 
corresponding object in the conformal turbulence, we need to compare the dissipation calculated 
by the viscosity term and nonlinear term in the Navier-Stokes equation. The correlation functions 
which we need to use is, 

<Uj{z)52m> = L-2^x|^|-2A^-2A,+2A,-6_ (52) 

If we denote the short distance cutoff by 5, the dissipation rate of the enstrophy can be either 
computed as, 

by using the viscosity term, or can be alternatively given by, 

by using the nonlinear term. In order to balance those two, we have to impose, 

^^52A^-4A,+2A^^ (55) 

The similar expression when we have vanishing one-point function is discussed by Hanany et. al., 

^^^2A^-2A,_ (5g) 

By employing these relations, one may obtain the dependence of r/ on the viscosity as 

^ 1-7 2A^ - + 3 

' " = "^ = A,-2A, + A, - ^''^ 

If we use the vanishing one-point function condition, 71 should be replaced by, 

2A^ + 3 

7 = 72 = -r (58) 

In either case, if we employ the Kraichnan-Leith-Batchelor type condition, one has to impose, 

7, = 1, i = l,2, (59) 
and if we use Saffman type condition, we get, 

li = \. 1 = 1,2. (60) 

On the other hand, the parameter /3 defined in ( p!9| ) may be expressed by the conformal 
dimensions as 



Q!i + 5' 

/?i = (2A^ - 4A<^ + 2Ax)-\ 

h = (2A^-2A<^)-\ (61) 
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by comparing with (55-5^). In this way, one can relate the dimensions of CFT primary fields 
with all the parameters which appeared in the previous section. 

4 CFT models 

In order to compute concrete values for the critical exponents, we need to solve (generalized) 
Polyakov's conditions (|57| ) or (^) explicitly. This program is rather hard to be carried out 
analytically, hence, we have only solutions obtained from numerical computation. Previously, 



some of the CFT models of KLB type was found by |22| for vanishing one-point function case 
and by Q for non-vanishing one point function case. In those papers, although CFT models are 
derived, the characteristic features which are common for those models are not discussed. In this 
paper, we study more extensive study for the generalized conditions (^7|) or (^8|) to get insights on 
the generic features of models which meet these conditions. We surveyed generalized Diofantian 
equations ( |59| ) (|6C|) (we shall call them as Pi, and Si (with i = 1,2) in the following) up to 
p < 1000 and q < 100. As expected we got many CFT candidates, 1658 models for condition PI, 
78 models for condition P2, 67 models for condition SI and 75 models for condition S2. Later we 
prove that there are infinite number of solutions for any conditions of type (^). 

Since we have several CFT models, we can not determine the exponents in the two dimensional 
turbulence uniquely. In figure la-d, we show the histogram of number of CFT models against the 
exponent of energy spectrum or each of the conditions. Where as in figure lb, Id, we have rather 
flat spectrum in the range —5 < a < —3, we have broader (—5 < a < 1) and localized spectrum 
(peak located around a = —4.5, —3.5, —2.5, —1.5, —0.5, 0.5) for the modified condition. We 
can also notice that a significant portion (82%) of the CFT models has their exponent in the 
range — 5 < a < —4. 

We remark that in the course of computation, we need to determine operator (j) as the lowest 
dimension operator in the OPE tp x ip. By this process, we can classify the obtained models into 
two groups. 

The first group, which we call category A below, is characterized by the fact that (p has already 
the dimension very close to the lowest dimension Aq of the system. In this case, the dimension 
of X is very close to (p since both of them are roughly same as Aq. The situation remains same 
even if we take higher OPEs. Most of the solutions which we obtained belong to this category. 

On the other hand, in the second class of solutions (category B), neither (p nor x have dimen- 
sion closer to Aq. In this case, in every step we take OPE, we meet lower dimensional operator. 
Since Category B is simpler and can be classified completely, we start with it first. 
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4.1 Category B 

The most interesting feature of solutions of this class is that we can find an infinite number of 
solutions for the condition of type ( |34[ ) for any fractional 7. As far as we have observed, the 
solutions of those type always take the following form, i.e. 

{p,q) = {pi{n-lf,q), 

il) = (n, 1), 

4> = (2n-l,l), 

X = (3n-2,l). (62) 



It is obvious from (35) that there appears only primary fields of type (r, 1) in the OPE of primary 
fields with index s = 1. On the other hand, the conformal dimension of fields (r, 1) decreases as 
r increases as long as r < pi/q{n — 1)^. One may confirm that (j) is the minimum dimension in 
OPE ip X ip and so is x in OPE ip x (f) if 



Sy + 2 + V9j/2 + Ay 
n > , y = q/pi. (63) 

Also, for sufficiently large n, one may easily convince oneself that neither (j) nor x have dimension 
closer to the minimum dimension of the system. 

Now if we put those combinations into (|3^), we get. 



1-7 _ 2A^-A0 + 3 6pi-q , , 

^ ~ 1^ 7 = n — — — 7— = • (64) 

The right hand side does not depend on n. For any fractional value for 7, we can adjust pi and 
q such that (]6^ ) can be satisfied. For example, to achieve KLB type condition (7 = 1), one may 
put 

Pi = l, g = 3. (65) 
For Saffman type condition (7 = 1/3), one may get. 

Pi = 2, q = 9. (66) 

These give the solutions to the Diofantian equation for n > 10 in KLB case and for n > 14 for 
Saffman case. This simple observation shows that there are infinite number of solutions for the 



condition (34) for any 7. 

One of the interesting feature of these solutions is that all of these models have same energy 
spectrum for each 7. Indeed, one may derive another cancellation of n, 

a = 4A^-2A<A + l = -^^^. (67) 

Pi 
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One can remove pi and q from relations (^) and (| 



a = (68) 
7 + 1 

The parameter /?, on the other hand, can be similarly determined by, 

^=^ = ^- (69) 
a + 5 6 

For KLB, one has (a,/3) = (—2,1/3) and for Saffman type condition, one has (a,/?) = 
(—7/2,2/9). The latter model seems to fit rather well to the observed value for the energy 
spectrum —4 < a < —3. 

4.2 Category A 

Although CFT models in category B has rather interesting properties, the majority (99%) of the 
solutions we obtained belong to this category. Although we can not solve the Diofantian equation 
in general, we can make several observations for this class. 

Interestingly, the various exponents of the models in this class is (although approximately) 
related to the central charge of the system. The dependence itself is different for the conditions 
that we impose on the conformal models. If there are some methods where we can independently 
measure the central charge, one may judge which condition is the appropriate one. 

In order to derive such formulae, it is essential to observe that 

1. The operator product expansion of ip field gives rise to a operator (p whose dimension is 
very close to the lowest dimension (Aq) of the system. 

2. The central charge of the nonunitary minimal model is linearly related to Aq. 

To demonstrate the first fact, we make histograms of the difference of the conformal dimen- 
sions between field (p and Aq (figure 2). We can see that deviation is less than 0.02. The second 
fact can be easily checked by recalling that c = 1 — ^^^^^^ and Aq = ^ ^Ipcf* ■ When p or q is 
sufficiently large, the first term in the denominator of Aq may be neglected and one finds, 

A„ = '-i^. (70) 



If we use the generalized condition (|57D , we get 



(l-7i)Ao-3 

A^- . (71) 



The energy exponent may be then given by. 



TP(i\ uoci -271A0-71 - 10 

E{k) = k\ ai = . (72) 

2-71 
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In particular, 

ai = -ll-2Ao (KLB), -8 - Aq (Saffman). (73) 
Similarly, the viscous scale index /3 in ([T^) can be obtained as 



for the condition (p^). From the KLB and Saffman's condition, /3 reduces to 

^ = -^u7r-T\ (^L^)' -lu7rT^\ (Saffman). (75) 
2(Ao + 3) 6(Ao + 3) 

Together with (^), it is clear that (|72|) gives the linear relation for the energy exponent and 
the central charge. We plot those relation in figure 3, to show that these are fulfilled rather nicely. 

4.2.1 Families of infinite solutions in category A 

Although the solutions of category A seem to occur randomly, there are several models that can 
be classified into several families of solutions. We would like to explain it by using the examples 
of condition PI. 

First of all, as Chung et. al. has already pointed out, there are many occurrences of solutions 
(in our case, 1347 cases out of 1658), with special combination of primary fields, 

V^ = (26,2), (A =(15,1), x = (30,2). (76) 

We would like to prove that this combination becomes the solution of conformal turbulence 
whenever 

29 31 

y g < P < yg, (77) 

with p and q coprime. Roughly speaking, this shows that there are q solutions of this type for 
each q. This fact explains quite systematic occurrence of solutions of this type. 



The proof is quite elementary. We must first observe (76) solves ( p6| ) for arbitrary p and q. 
This requirement is equivalent to the algebraic equation 

rl+rl-rl-l = 0, 
si + si -si = 0, 
nsi + r2S2 - rsss -7 = 0, (78) 

if V', (t>, X are labeled (ri,si), (r2,S2), and (r3,S3). 



The next step is to study when (p described in (76) becomes the lowest dimension operator 



in OPE. By explicitly comparing the dimension of (p with its neighbors, (13, 1), (17, 1), (45,3), 
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(45,5), (47,5), one finds that ( [77| ) is the desired condition. Under this condition, one can easily 
show that X = (30, 2) is indeed the lowest dimension operator in OPE ip x (p. (QED) 

For this family of solutions, since the ratio of p and q is already given (|77|), one can determine 
the range of the exponents. For example, the range of a is immediately determined by using 
a = 3x + ^ - 87, 

- 4.672 <a< -4.177. (79) 

Since p can take arbitrary (of course, with coprime condition) integer value in (|77|), the distribu- 
tion in this range is quite flat. This fact is quite consistent with figure la. 

There are two other families of solutions with similar property, i.e. they satisfy (|78|). The 
first one is, 

^^ = (55,4), <^= (101,7), x = (115,8), (80) 

and the second one is, 

■0 = (362, 26), (/>= (209,15), x = (418,30). (81) 

There are 54 occurrence of the first family and 2 occurrence of the second one. One can similarly 
show that they contain infinite solutions as their member. The expected range of a is given by 

- 4.765 <a< -4.740, 

-4.998 <a< -4.999, (82) 

respectively. 

Although 1418 of 1658 solutions are explainable in this way, there are still 240 solutions that 
we can not describe their origin. However, we may give some 'phenomenological' reasoning why 
there appear peaks around —3.5, —2.5, -1.5,-0.5, and 0.5. As it turns out, those peaks comes 
from the models where (p happen to be (17,1), (19,1), (21,1), (23,1) and (25, 1).^ In order 
that those operators are roughly the lowest dimension fields, the ratio of p/q = x should be 
17,19,21,23,25 respectively. Putting these values in (^), we get, 

a = -3.47, -2.47, -1.48, -0.48, -0.52, (83) 

respectively. Of these classification, there are 180 models. 

To summarize, although CFT gives infinite number of solutions for Polyakov's criterion, we 
can extract quite a lot information from them. Especially, we can determine the relations between 
the exponents and have several models which have their a in the special range. 
■^We remark that the biggest peak around —4.5 comes from models with — (15, 1). 



18 



5 Higher Order Statistics 



As we saw in the previous sections, the primary fields tp, </)■■■ are related to the observable in 
two dimensional turbulence through the energy exponent. Now some natural question arises. In 
general, there are many other primary fields in the minimal models. How can we extract some 
informations for those fields? Here we propose that one direction we may take is to go further 
to calculate higher correlation functions. As long as the assumption of the local equilibrium is 
correct, if one makes the dimensionless combination such as, < tp'^ > / < ip'^ >'^/^, it should be 
a dimensionless constant and there should be no dependence on the Reynolds number. In CFT, 
this is not true because there are singularities in the short distance expansion. 

In the correlation function, < tpip-'-ip >, the dominant contribution comes from the one 
point function of the lowest dimension in the OPE x [-0] • • • [ip] . This gives the L dependence 
of the correlation function. On the other hand, the singularity coming from the short distance 
expansion of tp fields can be expressed as power of the short distance cutoff 6. Let us denote (pn 
as the lowest dimension operator in the nth order product and A„ as the dimension of (pn- We 
may evaluate the n point function in terms of them, 

It follows that the non-trivial dependence on 5 for the dimensionless combination. 

Here /?„ is the exponent for nth order products: 

Pn = -P{2An-nA^). (86) 

The simplest choice F3 and -F4 corresponds to what is usually called "skewness" and "flatness" 
in the fluid dynamics community. We shall denote the exponents for n = 3, 4 as /?5 and Pp in 
the following discussions. 

5.1 Category A 

For this category, one may substitute A„ = Aq in (|8^). Using ([75|), one obtains 



(5k = -p3{2Ak-kA^) 
Ao(fe-2)(7-2) 



2(3 + Ao) • ^''^ 
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5.2 Category B 

Skewness and Flatness may be computed by using, A3 = A-^ and A4 = A(4„_3 1), which is vahd 
for n > 2y + 1 + \J ^y^ + y. Amazingly, even the dependence on parameter 7 disappears. 

Us = -/3(2A^ - 3A<^) = -| 

/3f = -/3(2A4 - 4A<^) = -4. (88) 

In general, if the minimum dimension operator in /cth product of 'i/' is given by (fcn — A: + 1, 1), 
the A;th exponent is 

^(2 -A;) , , 

A = (89) 

In this way, we find that the solutions for category B gives unique prediction for the higher 
exponents. 

6 Casimir Invariants 

Up to now, CFT is supposed to describe the time independent aspect of the two dimensional tur- 
bulence. However, the use of Navier-Stokes equation gives a simple explanation time dependence 
for the decaying turbulence^. In particular, we would like to discuss the decay rate of Casimir 
invariants J7„ in detail to compare the results from CFT to our numerical computation. One may 
easily find it as, 

As usual, A„ is the minimum dimension in the OPE of n ^ fields. On the other hand, a 
dimensional argument shows 



.. = ^l^^y , (91) 



U i± / L\Oln 

where a„ denotes an anomalous dimension for r/„. Comparing ( pO| ) and (pi]), we find 

an = 2n(A^ + 1) - 2A„. (92) 
Using (|10D, (lit) and (H), (0) becomes 

r,^ = j^f{-/3(A^+l)+l}+l-/3(A„+l)|^:^^2/J(A^+l)"-2/3(An-l)-l_ ^^^^ 



^ This process seems to be similar to the deformation of conformal field theory by coupling the relevant operators 
to conformal invariant system discussed by Zamolodchikov Isif . In his discussion, the central charge c behaves as 
the Morse function in the theory space of CFT and increases its value in the process of renormalization group. In 
order to apply his argument to turbulence, we need to extend his argument to include the non-unitary theory. We 
would like to discuss this aspect in future issues. 
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Substituting (p^) into (p3[), we can estimate the decay rate of r2„ as 

r/nOc(-) t^", 

where the powers p„ and can be obtained as 

2/3(A^ + 1) + 1 - 7 2/?(l-A„)-7 
= ^ nH , 

/3(-2A^ + l) + l -2/3(1 -A„) + 2 

Qn = n -\ . 

^ 7-2 7-2 

On the other hand, the random phase approximation described in section 2 leads to 

p„ = (l-7)(l + ^), 
7-2 

n 

Qn = -1 + 



7-2 

We can see that this difference originates in the anomalous dimension. 
6.1 Category A 

Using (p5[), ([TII) and ([7^), the exponents p„, qn of decay rate can be obtained as 



2 , 2{-7i(Ao + l) + Ao-l} 

^" " Ao + 3''+ (Ao + 3)(7i-2) ^ 

Ao + 1 , 7i(Ao-l)+8 
— -n + 



They become 



Pn 
Qn 



Ao + 3 (Ao + 3)(7i-2) 

2(n - 2) 
~ Ao + 3 ' 
(Ao + l)n + Ao + 7 



for the KLB condition and 



Pn 

qn 



Ao + 3 



lOra + 4(Ao - 2) 
Ao + 3 ' 
5(Ao + l)n + Ao + 23 
Ao + 3 



for the Saffman condition. 
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6.2 Category B 

The decay rate for the higher Casimir mvariant can be found from (39). After some computation, 
we get 



2(7 + 1) 6(7 + l)(7-2) 3(7-2)' 

1 9 27^ + 77 - 4 4 + 7 , , 

n^+ ^/ — + ^ , ' . (101) 



2(7 + 1) 6(7 + l)(7-2) 3(7-2)- 
Unlike category A, these exponents grows as n^. They become 



for the KLB condition and, 



1 2 7 2 
Pr, = —n n-\ — , 

Fn ^ 12 3' 

Qn = X'''"^''"^' 



3 2 133 2 

p = -n n , 

8 120 15 

3 2^ 13 13 
Qn = -gn +^^n--, (103) 



for the Saffman condition. 



7 Numerical Analysis 

In this section, some statistical properties predicted by CFT are compared with the numerical 
data obtained by the direct simulation. The pseudospectral method is used, which efficiently 
computes convolutions of the nonlinear terms in the physical space and is fully dealiased by the 
two-third law. The modes are 256x256, the time integration is performed by the second-order 
Runge-Kutta scheme up to t = 10 and the width for time step is 0.002. We examined five different 
cases of viscosity, v = 0.0001, 0.0003, 0.001, 0.003 and 0.01 and main results for the first four 
cases are shown. Decaying turbulence created by (^ without external forcing is investigated. 

Temporal evolution of energy E, enstrophy Q. and enstrophy dissipation rate r] is shown in 
Figures 4a-c. Initial increase of r] is observed until t ~ 5 for v = 0.0001 and the decay of and r] 
is prominent after t = 7. Thus, data in the time interval 7 < t < 10 are used to fit the asymptotic 
behavior of ij: 

r,^{i^/Ent + to)i, (104) 

where t in ( p7| ) is replaced by t + to in order to get a more precise power from the finite-time 
integration. The time offset to is estimated by the asymptotic expression oi ^I/tj w a{t + to), 
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which is independent of q and the data 7 < i < 10. Numerically obtained values of parameters 
in ( |104| ) are shown in Table 2. 

Although the numerical power of t varies from —4 to —2.2 according to it is nearly constant 
when 0.001 < v < 0.003, which is obtained by the proper resolution and closer to that of Saffman 
than that of Batchelor. (See the energy spectrum in Figure 4d) The power p is obtained as —0.171 
by the two cases v = 0.001, 0.003 and the energy averaged in the time interval 7 < t < 10, which 
is close to but slightly larger than Saffman's value —0.133. Comparing ( |104] ) with (p7[), relations 
between p, q and 7 are obtained as 

P = 2li^. (105) 

7-2 

4 — T 

q = ^. (106) 

7-2 

Although 7 varies from 0.333 to 1.333 from the numerical data of q, the plausible value of 
p = —0.171 leads to 7 = 0.614 or 0.557. From (|2^ ) we obtain r/ ~ i>^~'^ and the power lies 
between 0.386 to 0.443. The obtained power locates between the KLB value and the Saffman's 
value 2/3. 

The temporal evolution of the inviscid invariant Qn and its dissipation rate r/^ is examined 
and shown in Figures 5a,b for 1^ = 0.0003 and even n between 4 and 20. For odd n, fi^ oscillates 
irregularly around zero and no clear dependence on z/ or n is observed. The adopted numerical 
method does not conserve even if z^* — because of truncation in the Fourier space, although 
energy and enstrophy are conserved. In fact, oscillations appear in i7„ in the case of the smallest 
viscosity and large n. However, we analyze the data of ry„ since the deviation is not monotonic 
and the decay dominates the oscillation for the most cases. 

Similar to the ??, the following asymptotic form of 7]n for even n is assumed: 

rin'x{u/E)P"{t + toy". (107) 

Figure 3a shows the powers qn versus n for the four cases. If there is an n dependence of qn 
predicted by conformal field theory like 

qn = Qo + nqi, (108) 

we may estimate them by the data shown in Figure 7. Table 3 shows estimated values of qo and 
qi. Constants go a^^e within the values of ( |99D or ( |100| ) predicted by category A. However, for qi, 
numerical values are smaller than theoretical ones. 
The other parameter pn, denoted by 

Pn=Po + npi, (109) 
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is also estimated by plotting log Cn/ log(z^/-E') where C„ is r/„/(t + to)''"- We find negative values, 
between -0.422 and -0.265, of po and, between -0.0145 and -0.00525, of pi the latter is very 
small. Since the constants Cn include a constant proportional to n?, the degree of coincidence is 
unsatisfactory. 

For skewness and flatness, we should determine what kind of variables is to be measured. We 
choose 



and 



s=- (110) 

F=^44 (111) 



since they show rapid growth in the initial stage t « 2. The former is related with production of 
palinstrophy. In Figure 7a, b maxima of S and F versus at the initial increase are plotted. We 
observe that both S and F increases as u decreases, as is observed in many experiments of three- 
dimensional turbulence |3^. Intermittency of turbulence is responsible for the tendency that 
skewness and flatness increase as Reynolds number becomes large. We observe that S ~ Re^'^^''^ 
and F « Re^ '^'^^. On the other hand, the predicted feature for S and F is the opposite, i.e., the 
decreases as Reynolds number increases. Actually this prediction contradicts with observations so 
far. Since the behaviour of higher order statistics is strongly dependent on what is the measured 
physical quantity, there may exist observables faithful to the present prediction. 



8 Conclusions and Discussions 

It is shown that the extended theory of two-dimensional turbulence characterized by the energy 
exponent and the viscous scale index can clarify relations between classical and conformal field 
theories. New classes of CFT solutions for two-dimensional turbulence are obtained and its 
statistical predictions are explored. Many of the obtained exponents of the energy spectrum lie 
between or close to the reliable values —3 and —4. The dependence of Casimir invariants on 
the energy, enstrophy and viscosity is shown and used to predict its decay rate. Dependence 
of skewness and flatness on Reynolds number are also revealed and the predicted tendency is 
contrary to intermittency effects in three-dimensional turbulence. This may imply that the 
quantities to be measured are not the usual ones like doj/dx. 

There are some possibilities for further study of higher order quantities. For flatness, behaviors 
of other quantities, e.g., iIj, uj, can be investigated if the integration time is so large that isolated 
coherent vortices emerge. To perform more precise numerical investigation on < >, we can 
invoke the Zeitlin's finite-mode system, which was studied numerically in detail by Hattori 
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without viscosity. In this system the efficient FFT method is not available and its simulation 
requires much CPU time for integration. 
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10 Figure Captions 

Figure 1. Histogram of the CFT solutions versus energy exponent, by using KLB condition (a), 
and Saffman condition (b). In both of the case, we assume that the one point function is not 
vanishing. 

Figure 2. Histograms of the difference between the dimension of (j) and the lowest dimension of 
the system. We use KLB condition (a) and Saffman condition (b). Here we removed datum of 
category B. 

Figure 3. Relation between the central charge and and energy exponent. We use KLB condition 
(a) and Saffman condition (b). 

Figure 4. Temporal evolution of the energy (a), enstrophy (b) and enstrophy dissipation rate (c) 
and energy spectrum (d) a.t t = 10. Four cases with the viscosity v = 0.0001,0.0003,0.001,0.003 
are drawn by solid, dashed, dotted and dot-dashed lines. A dashed straight line in (d) shows the 
k~A energy spectrum. 

Figure 5. Temporal evolution of Vtn (a) and ??„ (b) for the case v = 0.0003. Values for n = 
4, 6, 8, • • • 20 are denoted in the order from below. 



Figure 6. Numerically obtained values of the powers qn in (108). Cases v = 0.0001, 0.0003, 0.001, 
0.003 are denoted by a circle, square, diamond and cross respectively. 

Figure 7. Numerically obtained maxima in the initial growth of skewness (a) and flatness (b) 
defined in (|Tl0|) and in ( pl| ) versus Re. 



11 Table Captions 

Table 1. First few solutions for CFT turbulence, (a) KLB condition (with non- vanishing one- 
point function) 
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(p>q) 


V' 


</- 


X 


a 


(33,2) 


Aio,i ~ -3.00 


Ai7,i 3.64 


Ai6,i ~ -3.64 


-3.73 


(43,3) 


Ai2,i 3.01 


Ai5,i ~ -3.09 


Ai4,i 3.10 


-4.84 


(44,3) 


Ai8,i ~ -2.99 


Ai5,i ~ -3.18 


Ai4^i ^ —3.18 


-4.61 


(46,3) 


A2o,i ~ -2.99 


Ai5,i ~ -3.35 


Ai6,i ~ -3.34 


-4.28 


(52,3) 


A25,i ~ -3.00 


Ai7,i 3.85 


Ai7,i 3.85 


-3.31 


(67,3) 


A36,i ~ -3.00 


A23,i ~ -5.09 


A22,i ~ -5.09 


-0.84 


(77,3) 


A43,i 3.00 


A25,i ~ -5.92 


A25,i ~ -5.92 


0.84 


(100,3) 


Aii,i ~ -4.10 


A2i,i 6.70 


Aai,! ~ -7.80 


-2.00 


(121,3) 


Ai2,i ~ -4.61 


A23,i ~ —7.73 


A34,i ~ -9.34 


-2.00 



(b) Saffman condition (with non-vanishing one-point function) 



(p,q) 


V' 




X 


a 


(97,4) 


Ai2,i ~ -4.03 


A23,i ~ —5.56 


A24,i ~ —5.57 


-3.99 


(74,5) 


Ai2,i ~ -3.08 


Ai5,i ~ —3.22 


Ai4ji ~ —3.21 


-4.91 


(86,5) 


A23,i ~ -3.33 


Ai7,i 3.81 


Ai7,i ~ -3.81 


-4.67 


(139,5) 


A43^i r 4.38 


A27,i ~ —6.45 


A27 1 —6.45 


-3.62 


(284,5) 


A95,i ~ -7.28 


A57,i ~ -13.70 


A57,i ~ -13.70 


-0.72 


(131,6) 


A32,i ~ -3.79 


A2i,i ~ -4.96 


A22,i ' 4.97 


-4.22 


(117,7) 


A22,i ~ -3.28 


Ai7,i 3.69 


Ai6,i ~ —3.69 


-4.72 


(299,8) 


Ai4,i 5.20 


A27,i 8.13 


A38,i ~ -8.85 


-3.52 


(196,9) 


A77,4 ^ —3.78 


Aio9,5 ' 4.96 


A87,4 ^ -4.96 


-4.22 



(c) KLB condition (with vanishing one-point function) 



(p,q) 


V' 


4> 


X 


a 


(21,2) 


A4,i ~ -1.14 


A7,i ~ -1.86 


Aio,i ~ -2.14 


-3.57 


(25,3) 


Aiij ~ -1.40 


Ag,! ~ -1.60 


Ag,i ~ -1.60 


-4.60 


(26,3) 


As.i 1.31 


Ag,i ~ -1.69 


Ag,i ~ -1.69 


-4.23 


(55,6) 


Ai4 1 ~ — 1.18 


Ag,i ~ -1.82 


Aio,i - -1.80 


-3.73 


(62,7) 


Ai3,i ~ -1.26 


Ag,i ~ -1.74 


Ag,l ~ -1.74 


-4.03 


(67,8) 


A28,3 ~ —1.38 


A25,3 - -1-62 


A42,5 ~ -1.62 


-4.51 


(71,9) 


A32,4 1-50 


A55,7 - -1.50 


A32,4 ~ —1.50 


-4.99 


(87,11) 


Ai6,2 ~ —1.51 


A23,3 ^ -1-49 


Ai6,2 —1.51 


-5.03 


(91,11) 


Ai4,2 ~ —1.40 


A25,3 ~ -1.60 


Ai6,2 - -1.59 


-4.61 



(d) Saffman condition (with vanishing one-point function) 
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(p>q) 







X 


a 


(17,2) 


Ae,! ~ -1.47 


Ag,! 1.65 


Ag,! 1.65 


-4.88 


(41,2) 


A4,i 1.32 


A7,i ~ -2.41 


Aio,i 3.29 


-4.27 


(46,3) 


A27 1 ~ —1.13 


Ai5^i ~ —3.35 


Ai5_i ~ —3.35 


-3.52 


(46,5) 


Al3,l r 1.43 


Ag,! 1.83 


Ag,i 1.83 


-4.74 


(53,6) 


Ai2,i ~ -1.45 


Ag,l r 1.74 


Ag,! 1.72 


-4.81 


(65,6) 


Ai7,i 1.35 


Aii,i 2.23 


All 1 ~ —2.23 


-4.42 


(91,8) 


A41 3 ~ —1.33 


A57,5 ~ -2.37 


A57,5 ~ -2.37 


-4.31 


(131,8) 


A62,3 1-08 


^49,3 ~ -3.61 


^82,5 ~ -3.61 


-3.31 


(79,10) 


A39,5 ~ -1-50 


A71 g ~ —1.51 


A39,5 ~ -1.50 


-4.99 



Table 2. Numerically obtained values of the time offset Iq, the power q in ( [104D and 7 in ( |106| ). 



V 




Q 


7 


0.0001 


24.3 


-4.03 


1.34 


0.0003 


16.8 


-3.61 


1.23 


0.001 


3.58 


-2.28 


0.438 


0.003 


1.21 


-2.36 


0.529 



Table 3. Numerically obtained values of the constants qo and qi in ( |108| ). 



V 


go (4 < n < 20) 


gi (4 < n < 20) 


go (12 < n < 20) 


q\ (12 < n < 20) 


0.0001 


-7.41 


-0.326 


-8.64 


-0.245 


0.0003 


-5.84 


-0.0997 


-3.43 


-0.245 


0.001 


-2.95 


-0.126 


-1.35 


-0.225 


0.003 


-1.86 


-0.451 


-1.05 


-0.503 
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